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1. INTRODUCTION 



That there are models in which symmetries, spontaneously broken at low (zero) tem- 
perature, are restored at high (finite) temperature was conjectured [1], and explicitly con- 
firmed [2,3] a number of years ago. Quantitative investigations of the Abelian Higgs model 
have exposed what has become known as the Unitary Gauge Puzzle [2-6] . This puzzle arose 
when attempts were made to estimate the transition or critical temperature for the restora- 
tion of the symmetry; in the unitary gauge, a different value for the critical temperature 
was obtained as compared to those obtained in certain other gauges. This difference has, 
in the past, been attributed to pathologies of the unitary gauge and the related question 
of the validity of the loop expansion. The purpose of this paper is to provide a detailed 
explanation for the discrepancy as well as a proof that identical results for the critical 
temperature can in principle be obtained in the unitary gauge and other gauges provided 
that one uses the correct definition of the critical temperature and that it is evaluated in 
a self-consistent expansion. 

In order to set the stage for the detailed presentation of this puzzle in the next section, 
we first briefiy review the usual approach to calculating the critical temperature using the 
effective potential [2]. This function, being defined for constant field configurations, enables 
the survey of a large class of translationally invariant expectation values of the field. The 
ground state of the system is identified with the global minimum of the effective potential, 
and spontaneous symmetry breaking is signalled by a non-zero vacuum expectation value 
(VEV) of the field. Under the assumption that the transition to the symmetric phase 
takes place in a continuous {i.e., second-order) fashion, the relevant VEVs must go to zero 
smoothly as the temperature is raised to its critical value. Thus, the "origin" must be a 
local maximum of the effective potential in the broken phase, and the condition signaling 
the onset of the second-order phase transition is that the origin be both an extremum and 
an inflection point of the effective potential. 

For concreteness, we consider a complex scalar fleld parameterized in terms of a real 
"Higgs" mode (pn and a real "Goldstone" boson (pa- The order-parameter, whose vanishing 
signifles the onset of the phase transition, is provisionally identified with the effective 
thermal mass-squared of the Higgs field: 
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The second derivative is evaluated at the value of the VEVs which gives rise to the mini- 
mum of the effective potential. Thus, as the temperature increases from below the critical 
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temperature, both rnjj^^^ and the VEV decrease towards zero. Above the critical temper- 
ature, rrijj^i^ is positive and the VEV is zero. The critical temperature is thus defined to 
be that temperature which causes 111%^^^ to vanish at zero VEV: 



This definition of the critical temperature, however, is problematic for gauge theo- 
ries, since the effective potential, calculated from off-shell Green functions at zero external 
four-momentum, is gauge dependent [2,7] and parameterization dependent [8]. Several 
approaches may be taken to deal with these dependences. One is to consider only objects 
which are manifestly gauge independent (in this context, the Vilkovisky-DeWitt effective 
action [9] is particularly suitable) , as the critical temperature defined by such means must 
be gauge and parameterization invariant [10]. Alternatively, one can look for signals for 
symmetry breaking in terms of expectation values of gauge invariant composite opera- 
tors [11]. Another approach is to try to extract gauge invariant information from gauge 
dependent objects like the usual effective potential [7]. In this regard the "Ward-like" 
Nielsen identities are useful in determining the gauge and parameterization dependence of 
quantities. For example, the Nielson identities have been exploited to show that symmetry 
breaking and restoration are gauge invariant phenomena [12]. Recently, they have been 
generalized to obtain a set of algebraic identities which describe the gauge dependence of 
all n-point functions derivable from the effective action [13]. As we shall demonstrate, 
these identities guarantee that the same result for the critical temperature in any gauge 
or parameterization can be obtained order-by-order in a self-consistent expansion. 

This paper is organized as follows. In Section 2, the unitary gauge puzzle shall be 
briefiy reviewed and previous analyses of it discussed. We shall emphasize results found 
by Ueda [4] and by Arnold, Braaten and Vokos (ABV) [6]. Section 3 discusses the gauge 
dependence identities and how they provide a formal resolution to the unitary gauge puzzle 
to all loop-orders. A two-loop calculation by ABV is then shown to corroborate these 
formal arguments. Further support for this resolution is presented in Section 4, where we 
consider a complex theory which also suffers from a "unitary gauge puzzle" for different 
parameterizations of the complex field [8] . Using the definition of the critical temperature 
discussed in Section 3, we demonstrate that the correct result is obtained in a large class of 
parameterizations. This class includes a parameterization corresponding to the one which 
gave rise to the original unitary gauge puzzle. Section 5 contains a brief restatement of 
our conclusions. 
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2. THE UNITARY GAUGE PUZZLE 

Consider a complex scalar field $ with a renormalizable self-interaction, m^|$p + 
(A/6)|$|^, respecting the inherent U{1) global symmetry of the system. By gauging this 
symmetry, a model of scalar electrodynamics is obtained. The model with = —fi^ is 
the standard Abelian Higgs model with spontaneous symmetry breaking: 

£ = {V,^)\V^^) + - _ I^^^^P. + jc,a^^^^ . (2.1) 

The gauge-covariant derivative couples the gauge field to the scalar with strength e: = 
[d^ + zeA^]$, and the ^7(1) field strength tensor is F^^, = d^A^ — d^^A^. The gauge fixing 
term has been left unspecified at this point. Before one is able to progress further, one 
must choose a parameterization for the complex field. There are two distinct options which 
shall be discussed in turn. 

2.1) Cartesian parameterization 

The complex field $ may be decomposed into its real and imaginary parts: 

^=^[^+i^. (2.2) 

Spontaneous symmetry breaking produces an asymmetric vacuum state about which the 
true quantum fields {^f,^f} fluctuate. Without loss of generality, the entire zero temper- 
ature VEV of the scalar field {i.e., po/\/2) may be assigned to the real component: 

% = Po+'^, '0 = 0+^/. (2.3) 

The vector field A^ fluctuates about a zero VEV. Note that the functional Jacobian cor- 
responding to the change of fleld variables from {$, ^*} ^ ^f} is unity and thus the 
path integral measure is trivial. 

Feynman rules based upon the Lagrangian written in terms of these shifted flelds may 
be derived in a straightforward manner. In order to do this, a gauge flxing term must be 
specified, but the results are independent of this choice. The order-parameter as defined 
by the curvature of the effective potential may also be expressed as the self-energy of the 
Higgs field ^ evaluated at zero external four- momentum. At one- loop: 
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The term is the bare tree-level mass-squared which is combined with the divergent 
one-loop vacuum self-energy and renormalized in the usual way. Explicit calculation of 
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the finite temperature contribution leads to, at leading-order, and in the high temperature 
limit, 

mk. = 2/^' + OiP-') ■ (2.5) 

The terms which arise from the scalar sector have coupling strength A, while those from 
the gauge-field loops couple with e^. It is assumed that A and are of the same order. 
Setting the expression for the order-parameter (2.5) to zero results in an equation which 
leads to the (correct) estimate of the critical temperature [2]: 



= 2//2 _ 



A 1 o' 



^ + 0(/3;'). (2.6) 



2.2) Polar parameterization 

An alternative approach to this analysis is to use the magnitude and phase represen- 
tation of the complex field: 

^=^pe'^/p°, (2.7) 
v2 

where ($) = po/\/2 is the VEV of the scalar field at zero temperature. Furthermore, all 
explicit dependence upon the Goldstone mode 9 is eliminated by choosing the "unitary" 
gauge: 

B^ = A^ + —d^e . (2.8) 
Spontaneous symmetry breaking affects the modulus field alone: 

p = Po+J. (2.9) 

Thus, only the "physical" degrees of freedom are present: a massive scalar Higgs boson 'f , 
and a massive vector meson B^. 

Feynman rules for the Lagrangian expressed in terms of these fields may then be de- 
rived. Note that the gauge for the vector field has been specified by (2.8). The model 
in this form has lost its property of being manifestly renormalizable due to the large mo- 
mentum behaviour of the vector meson propagator. The transformation of field variables 
gives rise to a non-trivial functional Jacobian in the path-integral measure, which may be 
incorporated into the Lagrangean in the usual manner. 



4host = -c(^l + ^Jc, (2.10) 

where {c, c} are scalar fields obeying Fermi-Dirac statistics. Notwithstanding these com- 
plications, the calculation of the one-loop self-energy of the Higgs field can be performed 
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in a straightforward manner. At zero external four-momentum, one finds in the high 
temperature hmit the following equation for the order parameter: 

= V - - \e'j, + 0{r') . (2.11) 

Setting this expression equal to zero produces an equation for the critical temperature, 

^^+0{p-'). (2.12) 

which does not agree with the equation (2.6) obtained in the Cartesian parameterization. 
This disagreement constitutes the unitary gauge puzzle. It is noteworthy that the contri- 
bution from the gauge field sector is the same in both instances. This suggests that the 
unitary gauge puzzle arises from the scalar sector of the Abelian Higgs model, rather than 
the gauge field sector. 

2.3) Previous Analyses of the Unitary Gauge Puzzle 

It has been suggested that the apparent non-renormalizability of the Abelian Higgs 
model in the unitary gauge and the concommitant question of the validity of the loop- 
expansion is to blame for the unitary gauge puzzle. However, Ueda [4] observed that 
by using an improved order-parameter [the second-derivative of the effective potential 
plus "H(P^ = — m^) — H(P^ = 0)"] it was possible to obtain the same result for the 
transition temperature in the unitary gauge as in other gauges. As we shall see in the 
next section, this corresponds to the pole-definition of the Higgs mass. In the Cartesian 
parameterization, the one-loop self-energy is independent of the external four-momentum 
and thus this additional term vanishes. However, in the unitary gauge the momentum- 
dependent vertices give rise to a non-zero contribution which exactly accounts for the 
difference between (2.6) and (2.12). 

More recently, ABV re-analysed the unitary gauge puzzle and determined two general 
criteria which must be met for a calculation of the critical temperature to yield a result 
which is meaningful. The first criterion is that the order-parameter itself must be the 
physical, gauge-invariant mass-squared of the Higgs field as defined by the location of the 
propagator pole. The second criterion is that the perturbative expansion and the ultimate 
evaluation of physical quantities under investigation must be self-consistently undertaken. 
They performed a calculation of the critical temperature in the unitary gauge in powers 
of T/Tc valid for low temperatures, by consistently extracting the leading temperature 
dependence in a loop-expansion. In this way, they argued that it was possible to ob- 
tain agreement between the Cartesian gauge and the unitary gauge values of the critical 



= 2//^- 



A 1 5 
12 + r 



6 



temperature, provided that certain unsuppressed terms of the form (T/Tg)^" for n > 2 
which arose in the unitary gauge calculation appeared with vanishing coefficients. Through 
careful analysis they were able to show explicitly that the O(T^) terms did indeed cancel 
among themselves up to two- loop-order. They were, however, unable to prove or verify 
that this property persisted for all higher-orders. Moreover, in the ABV method of analy- 
sis, it seemed that cancellations were required between one- and two-loop contributions to 
the self-energy. 

In the next section, a general framework justifying Ueda's approach shall be presented. 
Furthermore, in this context, the "miraculous" cancellations found by ABV demonstrate 
the consistency of the loop-expansion in the unitary gauge to two-loop-order. 

3. GAUGE DEPENDENCE IDENTITIES 

The off-shell effective action and the effective potential depend explicitly on the gauge 
and the field parameterization in which they are calculated. However, this dependence is 
controlled by a set of generalized Ward identities that can be derived from the definition 
of the effective action, or its Legendre transform, the generating functional for connected 
diagrams. A complete discussion of these identities, and further references, can be found 
in [13]. The generalized Ward identities yield an equation that determines the gauge 
dependence of the two-point function D^^ in a field theoretic model. In the absence of 
sources, the equation reads: 

^Dr^' = D'k'AX^ + DT^'AX^ , (3.1) 

where we employ a condensed notation in which {i, j, k, . . .} refer to all discrete field indices 
as well as to continuous spacetime coordinates. In the above, AXj is a quantity that can 
be calculated in perturbation theory. At this stage, however, (3.1) is exact as it stands. 
The crucial feature to note is that the change in the two-point function is proportional to 
the inverse-propagator itself. 

For the scalar mode in the spontaneously broken Abelian Higgs model at finite tem- 
perature, the identities acquire a particularly simple form in momentum space: 

A(/ _ ^2 _ n(p2. ^2^ ^ ^^2 _ ^2 _ n(p2. ^2^ y)^^^ _ (g ^) 

In the above equation, = 2|U^ is the renormalized tree- level mass of the Higgs field and 
a Minkowski metric is used. The self-energy H can in principle depend on both p, the 
external four-momentum, and on p, the VEV of the scalar field. Clearly, (3.2) implies that 
the propagator pole-position as determined by 

p^=m^ + U{p^-p^,T) (3.3) 



7 



is both gauge and parameterization invariant, provided of course that the corresponding 
AX is well behaved there. 

It must be emphasized that (3.2) is exact. Thus, the temperature dependent mass 
as defined by the propagator pole, and hence the critical temperature, can in principle 
be computed in either unitary or renormalizable gauges to yield the same result. This 
is not the case for the curvature of the effective potential, i.e., the inverse-propagator at 
zero external four-momentum, which is gauge dependent, in general. An exception to this 
arises when the self-energy is independent of momentum, as happens in the renormalizable 
gauges to one- loop-order. However, when the self-energy does depend on momentum, it is 
necessary to use the pole-position definition of the temperature dependent mass in order 
to obtain the correct result for the critical temperature. 

There is a line of argument that one might choose to follow in the hope of expediting 
the analysis which, unfortunately, leads to erroneous conclusions. At the precise value of 
the critical temperature, the physical situation demands that the VEV vanish and that 
the propagator pole occur at = 0. Thus, one is tempted to specialize (3.3), and solve 

= + n(0; 0, Tc) , (3.4) 

for the critical temperature. If (3.4) were the case, then it would always be possible to 
obtain the correct, gauge invariant, critical temperature from the effective potential. How- 
ever, this argument presumes an exact knowledge of the self-energy, whereas in practice 
it is necessary to evaluate it perturbatively. We shall now show that the critical tempera- 
ture obtained from the pole-position definition of the mass is gauge and parameterization 
invariant order- by-order in any self-consistent expansion (including the loop-expansion). 
This is not, however, true for calculations based on (3.4). 

There are two distinct and yet equivalent ways in which the dependence on the VEV 
can be consistently accounted for in the self-energy. The first is to treat the VEV as a 
parameter to be solved for perturbatively from the condition that the effective potential 
attain its minimum there. In this case, the self-energy is computed from the usual set of 
IPI diagrams. The second approach is to include "tadpole" diagrams in the calculation of 
the self-energy. In this case one demands that the tree-level tadpole vanish, thus obtaining 
the relation given below in (3.6). The parameter p appearing in (3.3) et.seq., is set to its 
tree-level value to enable consistent iterative solution for the pole-position. Higher-order 
tadpole loop contributions are responsible for the evolution of the VEV, and their inclusion 
in the self-energy takes these effects into account. In this section, we adopt the former 
approach because that is the framework employed by ABV, making it easier for us to 
utilize their results at two-loops in the unitary gauge (c./. Eq.(3.6) of Ref.[6]). However, 
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in the next section, the latter method is chosen, where only the self-energy needs to be 
calculated (with extra tadpole diagrams) rather than both the self-energy and the effective 
potential. 

One ordinarily employs a perturbative expansion that can be represented schemati- 
cally as follows: 



n = /n(i) + /2n(2) + ... , 

2 2 I 7 2 I 72 2 I 
P =P(0) +^ P(2) + -- - ' 



(3.5a) 
(3.56) 



Note that the self-energy starts at the first-order in the expansion, since we assume that 
the zeroth-order action is the classical action-^. Hence, the lowest-order scalar VEV must 
solve the classical field equations. For the Abelian Higgs model this leads to 



P(o) = 3m^/A = . 



(3.6) 



We also assume that renormalization (to absorb all divergences into the parameters of the 
theory) has been carried out consistently. In particular, the temperature dependent terms 
in the perturbative expansions considered here are all finite. 

In (3.5a, 6), / may be any parameter. It may be a combination of coupling constants, 
or a formal loop-counting parameter, such as %. It is vital that all of the contributions 
to any given order in the expansion parameter be properly taken into account. If this is 
done, then the pole-position, as determined by (3.3) is gauge fixing and parameterization 
independent order-by- order in I. Equation (3.3) must be solved with care. The only 
self-consistent solution takes the form of a Taylor expansion in the arbitrary parameter I: 



— 2 — 2 I 7 — 2 I 7 2 — 2 I 

P =P{0)+lP(l)+^ P{2) + --- > 

—2 —2 1 7 —2 I 72 —2 I 

P = PiO) + I P{1) + I P(2) + -- - ' 



where 



—2 2 
P{0) = TO 



p2^)=n(i)(p2;p2^T) 

P(2) = ^(2)(/;p^T) 



2 2 ) 

p =m 

P =P(0) 



2 2 

p =m 

2 —2 
P =P(0) 



, ^2 ^ 



2 2 

p =m 

2 —2 
P =P(0) 



2 2 ) 

p =m 

2 —2 
P =P(0) 



(3.7a) 
(3.76) 



(3.8a) 
(3.86) 

(3.8c) 



^ In field theories that require resummation due to infrared divergences this is not quite 
true, but this subtlety has no relevance to the present discussion. 
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and so on. The successive terms in the expansion of (3.76) are determined by analysis of the 
effective potential. The identities (3.2) guarantee algebraically that P(2)j etc., are each 
separately gauge and parameterization independent. Note that the correct evaluation of 
p^-^-^ has required setting the self-energy on its bare mass-shell = m^, as opposed top^ = 
0. In renormalizable gauges, since the self-energy is independent of the external momentum 
to this order, this distinction is moot. In unitary gauge, the distinction is crucial, and 
accounts for the discrepancy between the correct value for the critical temperature and 
that obtained from an expansion based on (3.4). Indeed, the gauge dependence identity for 
the two-point function, expanded to one-loop (first-order in the perturbative parameter) 
yields 

An(i)(0,p^o)>^) = -m'AX(i) , (3.9) 

which is not equal to zero in general. However, AX^^) = for renormalizable gauges and 
thus the method based on analysis of the effective potential works at one-loop in this class 
of gauges. In unitary gauge the self-energy is momentum dependent, and hence gauge 
and parameterization dependent as well, except on-shell. For cases of this type, the pole- 
position and associated self-consistent expansion outlined above must be used in order to 
get an invariant estimate for the particle mass and hence the critical temperature. 

In the Abelian Higgs model, it was first shown by Ueda [4] that, in the high temper- 
ature limit, 

n(i) = -2(A/2 + 3e')— - ^— (3.10a) 

in the unitary gauge. Thus the one-loop correction to the pole-position as obtained from 
(3.86), 



coincides precisely with the value obtained in renormalizable gauges. Fortunately, the 
elegant unitary gauge calculations of ABV are available to verify that the identities work 
at two-loops as well. The leading-order, high temperature contributions to the two-loop 
self-energy sum to 



while the shift in the scalar VEV is determined to be 



pJ,) = -^(A/2 + 3e2)^, (3.116) 
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and thus it is straightforward to verify that 

= 0(T2) . (3.12) 

The terms of order which contribute to the second-order correction in (3.12), appear 
with higher powers of the couphng constants, viz., {A^, e^, Ae^}, and so they are of higher- 
order than are the terms contributing to p'^-^-^ in (3.106) above. 

Note that throughout this section we treat I as a foop-counting parameter. As pointed 
out by ABV, in unitary gauge, the loop-expansion corresponds to an expansion in T^. In 
principle, therefore, the two-loop contribution could generate a large correction O(T^) to 
the one-loop result O(T^). The fact that this does not occur is a direct consequence of 
the gauge dependence identities, which guarantee that the same result to any order must 
be obtained in unitary gauge as in renormalizable gauges. In renormalizable gauges it is 
manifestly clear that the net O(T^) contribution must vanish at two-loop-order. Note that 
the expression (3.8c) for p^^) involves both 11(2) and derivatives of 11(1). ABV made this 
observation and interpreted it as a cancellation between one-loop and two-loop terms. They 
were led to conclude that generic calculations in unitary gauge would require contributions 
from all higher- loops in order to agree with results obtained in renormalizable gauges. In 
light of the gauge dependence identities, their analysis shows rather that one must keep 
all relevant contributions to the required order in the perturbative (in this case, the loop- 
counting) parameter. 

4. THE COMPLEX MODEL 

In order to see more clearly the issue of parameterization dependence in the pertur- 
bative computation of the critical temperature, we consider an interacting complex scalar 
field sufi'ering spontaneous symmetry breaking whose global U{1) symmetry is not gauged. 
The Lagrange density is 

c = (a^$)t(a^$) + /x2|$|2 _ A(|$|2)2 _ (4j) 

We introduce a one-parameter family of parameterizations in terms of two real fields 
{ ^} and a real parameter e: 

^=^[% + i^, (4.2) 

^ = (l-e) }f> + e\pcos{%/po) , (4.3) 
'(j> = {l-e)% + e\f>sm{%/po) . (4.4) 

Inspection of (4.2—4) reveals that the Cartesian parameterization of Section 2.1 is obtained 
for e = 0, while the polar parameterization corresponding to the unitary gauge in Section 
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2.2 results when e = 1. For < e < 1, an "interpolating" parameterization ensues. Both 
fields { }f), ^} are chosen to have the same canonical dimension as the scale associated 
with the angular field is identified here with the zero temperature VEV. This restriction 
shall be relaxed in Section 4.1. 

The Lagrangian (4.1) may be expressed in terms of { \f), ^} by substitution of (4.2 — 4). 
The zero temperature VEV of the scalar field may be chosen to be concentrated entirely 
in the ^ mode, 

V. = po+y, % = 0+J. (4.5) 

Substitution of (4.5) into the { \f), ^}-Lagrangian corresponds to the usual shift of the 
fields by their VEVs. The presence of the trigonometric functions in (4.3, 4) leads to the 
appearance of vertices to all orders in {]f, ^}. However, only those with up to four legs 
are required to construct the self-energies to one-loop-order. The Lagrangian is 



/: = c 



constant + linear ^quadratic 



-\- f A- C A- C A- C 

' ^quadratic ' ^cubic ' ^quartic ' ^hi; 



ligher-order 5 



(4.6) 



where 



constant -|- linear 

C 

^quadratic 

C 

*^cubic 

r 

^'quartic 



^22 ^4: 2 if 

2^ /^o - ^Po - ™2P0 7 , 

-\ y [d,d^ + mf] y - ^ y + (1 - 6)mi] y , 
-h^p' y' - ^^(1 - ^) y y' + ^^^^(^M ^)(^'^) ' 



-\\^<l-e){d,'f){d'^W 



4 Ap^ 



4^ 



1 



4 Ap^ 



r 



(4.7) 



and 



2 2,'^2 2 2i'^2 

m^ = -li + -po , ma = -p + -Po • 



(4.8) 



12 



Propagators: 



Cubic Vertices: 



fc2 



-^i2(gi;g2,g3)5(E9) 912 = ^(1 - e) - ^eg2 • gs 



-9icc{qi;q2,q3)S{^q) gicc = 



Quartic Vertices: 




\ ✓ 

X 

✓ \ 



>: 

✓ 



-Ai(gi,?2,?3,?4)(^(Eg) Ai = A 



-Ai2(gi,g2;g3,g4)5(E?) A12 = A 



-><2{qi,q2,q3,q4)S{J2q) A2 = A 



^i^-<^) + ^<^-^)qi-q2 



(1 - 



qi-qs + qi- q^ 

+92 •93 + 92-94 



91 ■ 92 + <?i ■ 93 

+9i ■94 + 92-93 
+92 - 94 + 93 - 94 



-A2cc(9i, 92; 93, 94)5(X) 9) A2CC = -^^(1 - ^) 



Table 1: Feynman rules for the interacting complex scalar model expressed in terms of the 
interpolating parameterization. The solid, dashed and dotted lines denote Higgs (type- 
one), Goldstone (type-two) and ghost field propagators, respectively. The Goldstone boson 
is massless as a consequence of the tree-level stability condition. For the vertices, momenta 
are positive incoming and are enumerated in an anticlockwise fashion beginning from 
the upper left corner. Each of the vertices contains an energy(frequency)-momentum 
conserving delta-function. 



The presence of the derivative interactions obscures the manifest renormalizability 
of the model. In addition, the field transformations from {$,$*} — > { y, ^} produce a 
non-trivial Jacobian term in the functional measure, which can be represented through 
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Grassman-valued ghost fields as 



4host = -c(l + -e]f- ^e(l -e)T + ---)c. 
\ Po Po J 



(4.9) 



The terms which have been omitted (represented by eUipses) are of higher-order in the 
fields and hence will not contribute to the one-loop self-energies. 

Feynman rules for this parameterization at finite temperature in the imaginary time 
formalism are collected in Table 1. Appropriate to this formalism, in this section we work 
with Euclidean four-momentum P^; the Minkowski mass-shell relation = vn? of the 
previous section translates here into = —w?. As discussed in Section 3, the viewpoint 
which shall be adopted here is to incorporate the effects of the evolution of the VEV by 
inclusion of tadpole graphs in the computation of the self-energies. By adhering to this 
prescription, only the self-energy needs to be computed. The parametric value of the VEV 
in (3.3) is set by the classical Lagrangian, see (4.10), in order that the perturbative series be 
well-formed. Discarding the linear term in the Lagrangian (4.6) corresponds to dispensing 
with the tree-level tadpole diagram, and leads to the tree-level stability condition c.f., 
(3.6): 

mi = ^ m\ = \pl = 2ij? ^ Po=^- (4.10) 

The collection of diagrams contributing to the self-energy of the Higgs mode ( ]f) is 
assembled in Figure 1, with the corresponding expressions found in Table 2. 




Figure 1: The self-energy of the Higgs field to one-loop-order as expressed in terms of the 
interpolating parameterization. 

The effective mass-squared of the Goldstone boson shall also be explicitly calculated 
so as to provide a consistency check on the application of the method and its results. The 
set of self-energy diagrams for the Goldstone mode (^) is found in Figure 2, and their 
expressions are listed in Table 3. 
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Diagram 



-n(ii)(p) 



Naive Order 



al 

a2 

bl 
b2 
bg 
cl 

c2 

eg 



Jdk -^^A{k) 
Jdkl^^ik) 

^(^)} 

fdk-Ae 



Jdk 
Jdk 



2 ^2 



+ lower order terms 



1 ^2 



rp2 

rpA _j_ 2^2 

rp4 _j_ 2^2 

2^4 _j_ 2^2 

rp4 



Table 2: The self-energy contributions to the (11) two-point function at one-loop order 
which arise in the interacting complex scalar model as expressed in the interpolating pa- 
rameterization. Here, J dk = J j0^- 



al) 




a2) 



ag) 



\ I 



bl) 



b2) 



bg) 



cl2) 



Figure 2. The self-energy to one-loop-order of the Goldstone mode, as expressed in terms 
of the interpolating parameterization. 

In Tables 2 and 3, the leading naive degrees of temperature dependence associated 
with each diagram are given. We note that the vacuum contributions which are quartically 
divergent (associated with the O(T^) pieces) all cancel. 

To evaluate each of these self-energies, the field is chosen to be at rest with respect to 
the heat bath. In the high temperature limit, the following results hold: 



dk 



1 

FT 



thermal 



J{(3m) 
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Diagram 



-n(22)(p) 



Naive Order 



al 

a2 
ag 
bl 

b2 

bg 
cl2 



i(l-e) + Ke2 



/d/c<|^6(l-6)-A 



2ml 



Jdk 



■e{l-e) 



e(l-e) 



Jdk 



fdk- 



3mf 

^e(l -e) 



2P^ 



3m' 



e(l-e) 



/ dkX^e^iP ■ kf "Aik) ^{P - k) 



rp4 _j_ rp2 

rj-\A _|_ rj-\2 

rj\2 
rp2 



Table 3: The self-energy contributions to the (22) two-point function at one-loop-order 
which arise in the interacting complex scalar model as expressed in the interpolating pa- 
rameterization. 



dk 



{P-kf 



dk 



{k"^ +w?){{P -kf + w?) 
(A;2 + m2) ((P - A;)2 + M2) 



thermal 
at rest 



thermal 
at rest 



where 



dx- 



X 



1 



(4.11) 



Vx^T^ eVa;^+c^ - 1 
Moreover, in the high temperature limit /?m — > 0, 



(4.12) 



(4.13) 



This limit, as well as PM 0, is already assumed in equation (4.11). 

The order-parameter for the symmetry restoring phase transition is the effective mass- 
squared of the Higgs field. 



at rest 
high temp, 
on-shell 



(4.14) 



Note that renormalization is assumed to have taken place, and mf represents the physical 
value of the zero temperature Higgs boson mass-squared. In the high temperature (weak 
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coupling) limit, we find 



m 



leff 



- 



A 



1 - 



1 + 



i2e - 1) 



(4.15) 



on-shell 



Thus, on-shell at = — mf the coefficient of the e-dependent term vanishes identically, 
and the final expression for the order-parameter is 

A 



m 



leff 



= 2n' 



+ 0{/3-') . 



(4.16) 



Equation (4.16) is valid for the entire e- family of parameterizations and yields the correct 
critical temperature (c./., (2.5 — 6)) in all cases. This is the central result of our analysis 
of the interacting complex scalar model and has provided a practical illustration of the 
gauge-dependence identities. 

Note that in the e = limit of the general result (4.15) the term in square brackets 
vanishes, and the correct physical one-loop result for the value of the order-parameter 
(4.16) then follows independently of the application of the on-shell condition. This is due 
to the absence of derivative interactions in the Cartesian parameterization. Thus, the 
method of analysis based upon the effective potential approach as described in Section 2 
works in the Cartesian case to one- loop-order. Also note that for the polar case e ~ 1, 
evaluating the order parameter off-shell at = results in the (erroneous) result of the 
unitary gauge calculated from the effective potential alone: 



leff 



2//2 - 



A 



12/32 



(4.17) 



P=(0,0) 

As a test of the validity of the method and the consistency of the calculations, consider 
the expression for the effective mass-squared of the Goldstone mode: 



™2 _ 2 , rT(22) ,p2N 

'''■2eff ~ ''''2 ' J-J^TOTALl^-f^ ) 
thermal 



at rest 
high temp, 
on-shcll 



(4.18) 



We expect that the Goldstone boson is massless in the broken phase below T^, and indeed 
collecting all terms in the high temperature limit yields 



m 



2eff 



0- 



A 



12/32 [ Ap^ 



re(5-e) 



(4.19) 



which vanishes at = —m^ = 0, independent of e. Thus, the Goldstone boson remains 
massless at one-loop to the order at which it is calculated. 
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4.1) The Case of an Arbitrary Scale For the Goldstone Mode 

The calculation described above was undertaken with a parameterization which car- 
ried an explicit dependence on the VEV at zero temperature. In this subsection, a two- 
parameter family of parameterizations will be constructed in which the scale p is not 
identified with the VEV of the field: 

(4.20) 

e) V + e Vcos(^/p) , (4.21) 
e) ^ + e V)sin( Vp) • (4-22) 

This can be written in terms of the VEV po by introduction of a dimensionless parameter: 

p = ?7Po , < ?7 < oo . (4.23) 

A frequently occurring combination of r] and e shall be denoted by 6, where 

h = r} + e-r}e = r}+{l-r})e = e+{l-e)r} . (4.24) 

Although not necessary in this context, it is possible to choose a parameterization for 
which the standard normalization of the kinetic term for the Goldstone field results; this 
can be accomplished by using 

^ = (1 - e) I/, + e ^ cos( %/poh) , (4.25) 
= i^Lljh % + e\t>sm{ %/pob) , (4.26) 

rather than the parameterization of (4.21 - 4.22). VEVs are then assigned to the fields as 

^ = po+y, (4.27) 

The Feynman rules for this parameterization are collected in Table 4. Fortunately, no new 
vertices are required beyond those listed in Table 1, but the vertex factors are more com- 
plicated. As before, the tree-level stability condition (4.10) arises again when elimination 
of the linear term in the Lagrangian is imposed. 

The same diagrams listed in Figure 1 contribute to the self-energy in this more general 

parameterization. Evaluation of these diagrams in the high-temperature limit leads to the 

following expression for the order-parameter: 

+ 0{/3-^) . (4.28) 



'0=(1- 
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Propagators: 



Cubic Vertices: 



fc2 



-^i2(gi;g2,g3)5(E9) 912 = + ^[gi • (52 + 53)] - ^52 • gs 

-«/icc(qi; 52, g3)^(E 9icc = 



Quartic Vertices: 




X 



-Al(9l,92,93,g4)(^(E?) Ai = A 



-Ai2(gi,g2;g3,g4)5(E9) A12 = A 



-A2(gi,g2,g3,g4)5(Eg) A2 = A 



Ib^-e _ 4£(l-6) 
1 2e(l-e) 



gi • gs + gi ■ g4 
+g2 • gs + g2 • g4 



2e 



J7;2g3 • g4 



-A2cc(gi,g2;g3,g4)(^(Eg) A2 



Ap2 62 
P 

Ap264 



e(l-6)(l+77) 



gi • g2 + gi • gs 
+gi • g4 + g2 ■ gs 
+g2 • g4 + g3 • g4 



Table 4: Feynman rules for the interacting complex scalar model expressed in terms of the 
modified interpolating parameterization. The notation and conventions are the same as 
those used in Table 1. 



There are several things to note about this equation. The first is that setting the external 
four-momentum on-shell, = — mf, causes the coefficient of the e- and 6-dependent terms 
to vanish. Thus, the correct physical result is obtained for all values of the parameters. 
The second feature of (4.28) is that in the limit in which the scale is chosen to equal 
the VEV, viz., 77 = 1 and hence 6=1, then the result reduces to the one obtained 
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earlier in (4.15). The third aspect of (4.28) we mention is that its structure guarantees 
that the correct physical result is obtained to this order both on- and off-shell in the 
Cartesian parameterization, e = 0, as expected. Finally, it is noted that when the polar 
parameterization, e = 1, is chosen, then b = 1, and the result is manifestly independent of 
the scale chosen for the angular Goldstone field. This fortuituous independence enabled 
the unitary gauge puzzle to be more precisely formulated when it was first discovered. 

Finally, the Goldstone boson's thermal self-energy corrections may be calculated from 
the diagrams of Figure 2. The result is: 



where G{e, ry, b) is a function of its parameteric arguments. Thus, as was the case in (4.19), 
evaluating (4.29) on-shell at = shows that the Goldstone boson remains massless at 
one-loop-order independent of the precise details of the parameterization. 



We have shown that a complete resolution to the so-called unitary gauge puzzle is 
obtained by examining the Ward identities that govern the gauge and parameterization 
dependence of the quantities under consideration. These show directly that the discrepancy 
between the unitary gauge and renormalizable gauge results for the critical temperature 
derive not from any intrinsic pathologies of the unitary gauge, but from the momentum 
dependence of the self-energy. We argued that gauge invariant results for the critical 
temperature can be obtained order-by-order in any self-consistent perturbative expansion 
in any gauge, including unitary gauge. This was verified explicitly at the two-loop level 
to order in scalar QED, and at the one- loop level in a family of parameterizations 
in complex 0* theory. The latter calculation also stresses the important, but sometimes 
forgotten, fact that the unitary gauge puzzle is concerned more with parameterization 
invariance than with gauge invariance, although the two issues are deeply related. 

It must be stressed that our analysis does more than re-interpret the two-loop calcu- 
lations of ABV. First of all, our arguments, admittedly formal, provide an all-order proof 
that the cancellation observed by ABV in unitary gauge at the two-loop order occurs to 
all orders in perturbation theory. Secondly, our analysis provides criteria that can be used 
to determine whether or not a given perturbative calculation should be expected to yield 
gauge invariant results. In particular, it is necessary to start at lowest-order from a gauge 
invariant quantity (such as the pole-position of a physical propagator), and perform a 
self-consistent expansion. In the present paper the loop-counting parameter was used, but 




on-shell 



+ 




(4.29) 



5. CONCLUSIONS 
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of course any convenient parameter may be used (coupling constant, temperature, etc.). 
It is, however, important to note that although any such self-consistent expansion will 
guarantee gauge independence order-by-order, the question of accuracy must be addressed 
independently. The weakness of unitary guage is that this question is obscured due to the 
presence of non-renormalizable interactions introduced by the non-linear reparameteriza- 
tion. The utility of the gauge dependence identities lies in the fact that they guarantee 
that the error estimate can be done in any convenient gauge, such as renormalizable gauge, 
while the actual calculations can be performed in any other gauge, because the same an- 
swer will be obtained order-by- order in all gauges. Thus, for example, the accuracy of the 
leading one-loop result for the critical temperature in renormalizable gauges implies that 
higher-loop contributions in unitary gauge will also be negligable. 
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